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ABSTRACT 


Two  color  time  of  flight  mass  spectroscopy  is  employed  to  study  the  van 
der  Wauls  (vdW)  clusters  of  benzenefN^}^  (n<8),  benzenetCOg)^  (r<7),  and 
benzene(CO)^  (n=l,2)  created  in  a  supersonic  molecular  jet.  Potential  energy 
calculations  of  cluster  geometries,  normal  coordinate  analysis  of  vdW  vibra¬ 
tional  modes,  and  calculations  of  the  internal  rotational  transitions  are 
employed  for  the  assignment  of  the  benzene(solvent)^  cluster  spectra, ^n  the  0|’ 
and  6^  regions  of  the  benzene  1B2u  *  !Alg  transition.  The  respective  vibronic 
and  rotational  selection  rules  for  these  clusters  are  determined  based  on  the 
appropriate  point  groups  and  molecular  symmetry  groups  of  the  clusters.  Good 
agreement  between  the  calculated  and  experimental  spectra  is  obtained  with 

ic  >'■ 

regard  to  the  vdW  vibrational  and  internal  rotational  modes.  / The  solvent 
molecules  rotate  nearly  freely  with  respect  to  benzene  about  the  benzene- 
solvent  bond  axis  in  the  benzene(solvent ) j  clusters  In  the  excited  state  a 
small  ca.  20  cm*1  barrier  to  rotation  is  encountered.  Studies  of  larger 
clusters  (n>2)  reveal  a  broad  red  shifted  single  origin  in  the  6^  spectra.  A 
linearly  increasing  cluster  energy  shift  is  observed  as  a  function  of  cluster 
size.  The  cluster  energy  shifts  are  not  saturated  by  one  solvent  molecule  on 
each  side  of  the  aromatic  ring;  several  solvent  molecules  effectively  interact 
with  the  solute  k  electronic  cloud.  Both  homogeneous  and  inhomogeneous 
nucleation  take  place  for  the  clusters  studied  depending  on  the  ratio  of  the 
solvent-solvent  binding  energy  to  the  cluster  binding  energy. 


I . 


INTRODUCTION. 


Supersonic  molecular  jet  spectroscopy  has  made  possible  studies  of  a 
wide  variety  of  weakly  bound  molecular  van  der  Waals  (vdW)  clusters  vdW 
clusters  only  slightly  perturb  the  properties  of  their  individual  constituents 
and  are  characterized  by  relatively  low  binding  energies,  large  i ntermo  1  ecu ! nr 
equilibrium  distances,  and  low  frequency  intermolerular  vibrational  vdW  modes 
A  number  of  cluster  systems,  including,  for  example.  vdW  complexes  of 
aromatic  molecules  (e.g.,  benzene  and  toluene)  wi-K  hydrocarbon  solvents. 
and  clusters  of  N-heterocyc 1 i c  solutes  {e.g.,  pyrazine  and  pyrimidine)  with 
both  alkane'*'1*  and  hydrogen  bonding  solvents  (e.g.,  water  and  ammonia),4  have 
been  studied  in  our  laboratory.  Cluster  geometries,  energetics  and  energy 
dynamics  have  been  explored  in  these  systems  yielding  insight  into  nucleation 
processes  and  solvation  geometries  in  both  gas  and  condensed  phase  systems. 

Our  recent  studies  of  the  vdW  cluster  torsional  mode  structure  in  benzene- 
methane,5  -deuteromethane,  and  -carbon  tetraf luoride6  clusters  reveal 
that  these  clusters  are  rigid  with  regard  to  internal  rotation  of  the 
cluster  constituents.  The  internal  torsional  motion  is  found  to  be 
oscillatory  and  constrained  by  an  orientat ional ly  dependent 
intermolecular  potential  with  a  barrier  height  of  the  order  of  the 
cluster  binding  energy.  Given  the  geometries  of  the  systems  studied, 
the  rigidity  of  the  clusters  can  be  understood  by  considering  the  nature 
of  the  solvent  rotational  dynamics;  a  low  barrier  one-dimensional 
rotation  of  the  solvent  molecule  around  the  three-fold  solute-solvent 
bond  axis  is  not  possible  due  to  the  orientation  of  the  solvent  molecule 
rotational  axes. 

In  this  paper  we  report  spectroscopic  results  for  the  benzene(N2) j  , 
benzene(CO) ^ ,  and  benzene(C02 ) j  clusters:  these  systems  exhibit  almost  free 
one-dimensional  internal  rotation  between  the  cluster  solute  and  solvent. 


Benzene  is  chosen  as  a  ehromophore  molecule  since  the  formation  of  small 
clusters  with  benzene  is  relatively  well  understood.  The  high  symmetry  of 
benzene  makes  the  interpretation  of  cluster  geometries  and  theoretical  calcu¬ 
lations  of  rotational  and  vibrational  modes  relatively  straightforward.  The 
experimental  spectroscopic  results  (two  color  time  of  flight  mass  spectroscopy 
(TOFMS))  are  interpreted  with  the  aid  of  calculat ional  modeling  of  selected 
cluster  characteristics,  such  as  geometries,  binding  energies,  and  vdW  vibra¬ 
tional  and  rotational  energy  level  structures.  Cluster  geometry  and  favorable 
orientation  of  the  solvent  rotational  axis  permit  one-dimensional  internal 
rotational  motion  for  these  clusters.  The  analysis  of  the  internal  rotational 
energy  levels  of  the  vdW  clusters  studied  bears  some  resemblance  to  that 

carried  out  for  substituent  groups  in  non-rigid  aromatic  molecules  such  as 

7  8 

toluenes  and  xylenes. 

The  second  part  of  this  paper  is  devoted  to  spectroscopic  studies  of 
relatively  large  benzene(N2)n  (n«8).  benzene(CO)n  ' n=l , 2 )  and  benzone(C09)n 
(n-£7)  clusters.  Studies  of  the  formation  (nuclear  ion)  of  large  clusters  are 
important  since  these  clusters  are  potential  model  systems  for  condensed 
phases  and  because  of  an  intense  interest  in  nucleation  and  growth  of  aerosol 
particles.  We  have  previously  determined1,2  that  small  (n=l,2)  vdW  clusters 
of  benzene  and  toluene  with  hydrocarbon  solvents  (methane,  ethane,  propane) 
are  formed  via  solute  collisions  with  solvent  dimers,  trimers,  and  larger 
aggregates.  Both  homogeneous  and  inhomogeneous  nucleation  take  place  for 
these  systems  depending  on  the  cluster  binding  energy  compared  to  the  solvent- 
solvent  binding  energy. 

Larger  vdW  clusters  with  n>3  have  not  previously  been  extensively 
studied  spectroscopically  due  to  the  difficulties  in  obtaining  their  resolved 
optical  spectra.  Consequently,  little  spectroscopic  information  concerning 


these  clusters  is  available.  On  the  other  hand,  other  techniques,  such  as 

eg.  electron  diffraction 1  10  ami  infrared  spectroscopy11  are  typically  used 

2 

for  studies  of  much  larger  aggregates  (ri>10  )  but  do  not  give  specific  infor¬ 
mation  about  the  cluster  energetics  and  geometries. 

II.  EXPERIMENTAL  PROCEDURES 

A.  Spectroscopic  Techniques. 

The  S |  -  S()  absorption  spectra  of  the  clusters  are  obtained  using  a 
pulsed  supersonic  molecular  jet  expansion  and  two  color  TOFMS.  The  details  of 
the  experimental  apparatus  and  techniques  employed  have  been  described  ear¬ 
lier.1,2  Two  synchronized  N’d'2  YAG  pumped  dye  lasers  are  used  to  probe  both 
the  origin  and  the  6^  regions  of  the  benzene  -  A ^ g  transition  in  the 

cluster.  Exciton  LDS  698  dye  is  used  for  the  pump  laser  with  the  output 

+  o 

frequency  doubled  and  mixed  with  the  1.004  urn  fundamental  of  the  Nd  /YAG 
laser.  Either  Exciton  LDS  698  (doubled  and  mixed)  ur  F  548  (doubled)  dyes  are 
used  for  the  ionization  laser.  The  energy  of  the  i  .uization  laser  is  typical¬ 
ly  lowered  to  the  ionization  threshold  to  limit  cluster  fragmentation  upon 
ionization.  The  spectra  are  recorded  using  a  boxcar  averager,  transient 
digitizer  and  computer. 

Typically,  35%  mixtures  of  the  solvent  gases  with  helium  as  a  carrier 
gas  are  passed  through  a  liquid  benzene  trap  at  room  temperature  and  subse¬ 
quently  expanded  into  the  apparatus  vacuum  chamber  by  using  a  pulsed  nozzle. 

A  nozzle  backing  pressure  of  100  psi  is  typically  used  and  the  chamber 
pressure  is  maintained  below  4xl0-6  torr. 

B.  Calculations  of  Cluster  Geometries. 

The  cluster  ground  state  geometries  and  binding  energies  are 

calculated  for  small  clusters  (n=l,2)  via  an  intermoieculur  potential  energy 

1  2 

minimization  procedure  employing  the  methods  previously  described.  ’  A 


Lennard-.Jones  (L-.J;  6-12-1)  potent  ia!  is  used  w  i  t.h  the  atom  atom  interaction 
parameters  previously  reported  for  similar  conciliations  The  potential  form 
and  calculation'll  algorithm  yield  consistent  results  with  regard  to  the  numbe 
of  cluster  configur.it  ions  observed  experimentally,  their  respective  binding 
energies,  and  qualitative  geometries  (symmetries). 

ITT.  THEORETICAL  ANALYSIS  OF  CLUSTER  VIBRATIONAL  AND  INTERNAL  ROTATIONAL 

MOTION. 

A.  Normal  Coordinate  Analysis  of  vdW  Vibrations  and  Vibronic 
Transitions . 

The  normal  coordinate  analysis  (N’CA)  of  the  cluster  vdW  vibrations 
is  performed  by  using  the  same  methods  as  described  in  our  previous  publica¬ 
tion'"'  on  vdW  cluster  vi bionic  structure.  The  GF  method  of  Wilson*3  is 
employed  in  which  the  cluster  constituent  force  fields  are  generated  by  using 
the  central  force  approximation  including  out-of-plane  motion  terms.  The 
intramolecular  force  constants  chosen  correspond  to  general  functional  group 
stretches  and  bends.  *  ^  The  i ntermo  1  ecu  lar  force  c<»:  stants  are  generated  as 
second  derivatives  of  the  intermolecular  L-J  potential  discussed  above.  The 
GF  matrix  is  numerically  diagonalized  and  the  3N-6  (3N-5  in  the  present 
instance)  nonzero  eigenvector  normal  modes  and  eigenvalue  energies  are  deter¬ 
mined  in  the  usual  fashion.  The  reader  is  referred  to  our  previous  publica¬ 
tion0  for  details  of  the  theoretical  procedures  employed. 

The  eigenvector  normal  modes  found  based  on  the  calculated  geometries 
of  the  rigid  clusters  (point  groups  c2v  and  Cs)  are  I.  the  vdW  stretch  <$z) 
which  entails  the  motion  of  the  solvent  molecule  in  the  direction  perpend i c  - 
ular  to  the  benzene  molecular  plane,  2.  the  vdW  bends  (h„  arid  b„)  with 
translational  motions  parallel  to  the  benzene  ring  (x  or  y  directions),  and  3 
the  vdW  torsions  (ty  and  tz)  which  transform  as  the  rotations  of  the  cluster 
constituents  about  the  y  and  z  axes,  respectively.  These  normal  modes  are 


shown  if  Fig.  1  for  the  .  i  s  <  ■  of  the  rigid  benzene (Ng)  j  cluster.  The  ■ 
torsional  morion  will  no.'  he  prt‘st>nt  if  the  solvent  molecule  can  undergo 
nearly  free  rotation  about  the  cluster  z  axis,  perpendicular  to  its  bond  axis 
instead  additional  j,.na!  levels  will  he  present  in  the  eigenstates  Of 

course,  the  f\,  rigid  iuster  point  group  is  no  longer  applicable  fur  ♦  he 
description  of  the  rotational  it-wls  and  an  appropriate  molecular  symmetry 
group  must  he  used  We  issume  that  ill  the  low-lying  i n t er mo  1  ecu  !  a r  eigen¬ 
states  other  than  t  are  neat  !y  h  muonic. 

The  appropriate  selection  rules  for  vihronic  transitions  in  the  rigid 
cluster  are  determined  by  analyzing  qualitatively  transition  moment,  matrix 
elements  using  the  adiahat.it;  approximation.  In  the  C 2V  symmetry  of  the 
benzene(N’/,)|  and  benzene) CO/, },  rigid  clusters,  the  vdW  vihronic  transition 
selection  rules  for  the  o"  region  in  the  absence  of  the  Herzherg-Tel ler  (HT) 
coupling  are  as  follows:  iv  =0,  il.  ±2.  for  tl  ■  totally  symmetric  vdW 

stretch  and  Av  -  0,  z2.  r4 .  ...  for  the  nontotal  ly  ..mmetric  vdW  bends  and 
torsions.  In  the  region  the  fundamental  of  the  nontotal  ly  symmetric 
torsion  {t7)  is  additionally  allowed.  This  torsion  is  vibronically  induced  in 
the  0°  region  if  HT  coupling  is  active.  For  the  Cs  symmetry  of  the 
benzenelCOlj  rigid  cluster,  the  vibronic  selection  rules  for  the  0°  region  are 
slightly  different:  =  0,  tl ,  ±2,  ...  for  the  vdW  stretch  S  bend  b  .  and 

torsion  ty,  and  iv  =  0,  ±2.  ±4,  ...  for  the  bend  by  and  torsion  t  ?.  In  rg 

symmetry  all  the  vdW  modes  are  capable  of  vibronic  coupling  and,  hence,  ^v  = 

0,  tl.  t2 .  ...  if  HT  coupling  is  present,  .Av  =  0,  tl,  t2,  ...  for  all  the  vdW 
modes  even  in  the  absence  of  HT  coupling  for  this  latter  cluster  in  the 
region 

B.  Internal  Rotational  Analysis  of  the  \.?  Torsional  Mode. 

We  assume  in  the  present  analysis  that  the  solvent  molecule  can 

rotate  in  onp- d i mens  ion  about  its  axis  of  largest  moment  of  inertia,  i.e  .  the 


axis  perpendicular  to  the  \'  (f  0)  mo leru  1  ar  bond  and  passing  through  the 
solvent  and  solute  molecular  center  of  mass.  The  rotational  axes  for  N,  and 
C0o  in  the  cluster  correspond  to  the  solute-solvent  bond  (cluster  z)  ixis 
(benzene  six  fold  symmetiy  ixis)  and  the  CO  rotational  axis  :s  found  to  form  a 
on.  15'  angle  with  the  benzene  six  fold  axis.  The  appropriate  mo  If  u!  n 
symmetry  group'"  must  be  used  :n  ‘tie  determination  of  the  rotational  mode 
structure  and  selection  rules. 

Analysis  of  the  interna!  rotational  transitions  in  these  clusters  is 

O 

performed  by  using  the  same  approach  as  described  in  our  previous  publication 
for  internal  rotations  in  non-rigid  aromatic  molecules.  The  energy  levels  and 
wave func t ions  for  the  mot ;  of  a  one  dimensional  hindered  rotor  are  obtained 

by  solving  the  Schrtxii tiger  equation 

•> 

[-B  _  -  V  ( 0 )  ]  ^{0)  =  Em  0/0)  ,  (1) 

50 

in  which  the  potential  form  for  the  benzene  clusters  considered  is  approxi¬ 
mated  by 

V(0)  =  -  V6  (1  -  cos  60) .  (2) 

In  eq.  (1)  B  is  the  rotational  constant,  (h 'Bit  cl),  I  is  the  moment  of  inertia, 
and  0  is  the  torsional  angle  for  the  N’9,  C09  and  CO  rotors.  The  solution  of 
the  Schrodingor  equation  for  a  free  rotor  (V(0)=O)  is  of  the  form 


(3) 


and  the  free  rotor  energy  levels  are 


E 


m 


m  B . 


(4) 


The  eigenvalues  and  eigenvectors  of  the  hindered  rigid  rotor  Hamiltonian  eq. 


(1)  are  solved  for  within  a  basis  set  of  21  free  rotor  eigenfunctions.  The 


parameter  is  adjusted  ! 
experimental  data. 


the  energy  levels  of  file  hindered  rotor  to  rhe 


The  energy  levels  1‘  ;r  m'  .it  inn  of  the  solvent  with  respei  t  to  the 


solute  .ire  labeled  ,ia  nrd  ,  n~  f  tie  m  quantum  numbers  of  a  one-d  i  mens  i  nnu  ! 


rotor.  The  symmetry  of  •  da  nil:  eiiari  is  given  by  the  appropriate  mole  alt- 


symmetry  group:*''  ,  f  the  benzcneO.'-,)  j  and  henzenelCf^ ) ,  cluster 


O 

(iij2  for  the  benzene ( CO ) ,  '  luster  {according  to  the  notation  in  ref  16).  The 


symmetry  allowed  rotational  transit  ions  within  the  0^  and  8  bands  are  deter¬ 


mined  via  the  nonvanishing  integral  rule  taking  into  account  the  correlation 


between  the  molecular  symmetry  group  of  benzene  and  the  appropriate  molecular 


symmetry  groups  of  the  lusters  The  results  of  this  analysis  for  different 


clusters  are  given  in  Table  r  The  a^  --  .tjff  { a  ^  --  a  ^ )  and  b2u  --  b,..  ( b  2 
--  b0)  transitions  are  riot  allowed  at  the  origin  but  are  allowed  at  fi*.  The 


2g  "  2g 


e0  --  Pol  ■  p<  „  ^  •  Pi.,  (Pi  --  'hi,  a 


In  ‘  1  '  •  ulg  2g  '  1 


a,  e o  )  and  e, 


b,u  (e1  --  bo)  transitions  are  allowed  in  both  spe<  • al  regions.  This  finding 


is  important  for  the  interpretation  and  assignments  of  the  spectra,  partic. 


lariy  in  i  he  f)^  region.  The  dependence  of  the  energy  of  each  internal  rota¬ 


tional  level  on  the  parameter  Vfi  is  illustrated  for  the  case  of  the 


benzene(No) ,  cluster  in  Fig.  2.  Similar  dependences  are  obtained  for  the 


benzene(CO'j  and  benzene(COp)-j  clusters.  In  the  latter  case  the  spacings 


between  individual  rotational  levels  are  relatively  smaller  due  to  a  larger 


moment  of  inertia  of  the  C0o- 


Some  of  the  allowed  transitions  may  not  be  observed  due  to  nuclear  spi: 


selection  rules  Moreover,  the  different  nuclear  spin  states  associated  with 


various  rotational  levels  result  in  hot  bands  which  cannot  he  removed  from  the 


spectra  by  normal  cooling  techniques.  For  the  benzene(N'., ) .  and  benzene  (CO , ! 


clusters  (Gp.  molecular  symmetry  group)  only  the  "g"  internal  rotational 
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levels  tire  populated  with  the  following  relative  statistical  weights: 


alg  : 

elu  ^2u  :  e2g  =  1  :  0  :  0  ;  2.  "g"  and  "u"  type  levels  are  not  >|  i  s 

a 

tinguished  in  the  G,2  molecular  symmetry  group  and  all  the  different  rota 
tinna!  symmetry  levels  of  benzene!  m)  ,  may  be  populated.  The  rel.t'jve  s'.it  is 
t  i  im  1  weights  .ire .  in  this  latter  vase,  -i,  •  e  ,  :  b,  :  e9  = 

12:1:2 

IV.  BENZENE  (SOLVENT),  CLUSTERS  SPECTRA.  CALCULATIONS  AND  ASSIGNMENTS 

A  similar  approach  to  spectral  analysis  consisting  nf  calculating  the 
internal  rotational  and  vibrational  transitions  is  employed  for  all  three 
benzene  (solvent),  clusters  studied.  Therefore,  only  the  assignments  for  the 
benzene (N?) ,  spectra  are  discussed  in  detail;  other  cluster  spectra  are 
analyzed  based  on  the  same  arguments. 

A.  Benzene (N2) , . 

The  results  of  calculations  of  the  group  state  geometry  and  vdW 
vibrations  of  the  benzene(V2),  cluster  are  shown  in  Fig.  1.  The  calculations 
yield  a  single  cluster  geometry  with  the  solvent  molecule  located  in  the  plane 
parallel  to  benzene  at  a  distance  3.3  A  above  the  center  of  the  benzene  ring. 
The  calculated  ground  state  binding  energy  is  -501  cm  .Vo  other  local 
potential  minima  are  found  in  the  calculations.  The  eigenvectors  of  the 
vibrational  vdW  modes  and  their  corresponding  energies  calculated  from  the  VCA 
ire  shown  in  Fig.  1.  Rotating  the  N2  molecule  by  an  angle  30*  around  the  z 
.is  (V2  rotational  axis)  results  in  an  increase  of  the  calculated  binding 
energy  by  only  ca.  1  cm"*.  This  indicates  no  barrier  to  rotation  in  the 
ground  state  and  suggests  that  the  N2  molecule  can  rotate  in  the  cluster  with 
little  barrier  to  rotation. 

The  spectra  of  the  benzene(N?) ,  cluster  in  the  and  6^  regions  of  ’he 
benzene  ‘R2u  “  '  A ,  ^  transition  are  shown  in  Fig.  3.  Such  spectra  are  typi 


ca  1 1  y  understood 5  based  on  allowed  origins  for  the  0*'  and  6^  regions  and  a 
of  assigned  vdW  normal  mo<ies  built  on  these  origins.  The  origin  in  the 
cluster  0°  spectrum  (if  allowed)  would  typically  exhibit  a  shift  in  energy 
from  the  benzene  0°(  tratisit  i .in  similar  to  that,  observed  in  the  0^  spectrum  of 
the  cluster.  For  the  benzene  i  N9 : .  *  luster  this  shift  equals  -0  cm  We  are 

net  able  to  understand  these  spe,  u  a  and  find  appropriate  assignments  based  on 
this  now  standard  approach,  however  The  spectrum,  at  least,  must  be 
analyzed  in  the  light  of  *  he  proper  group  theory  and  the  vdW  rotation 
vibration  energy  level  set  discussed  above. 

If  the  N2  molecule  in  the  benzene(.\’2)  j  cluster  were  freely  rotating  in 

both  SQ  and  Sj.  0alg  -  0a,g.  1 1  u  "  le!tr  2p2g  “  2H2g'  and  4e2g  ~  4p2g  1  z 
rotational  transitions  (replacing  the  t y  vdW  mode)  would  coincide  in  energy 
and  give  only  one  peak  at  the  origin.  The  first  two  of  these  transitions  are 
not  allowed  in  the  0^  region  due  either  to  vibronh  u*  nuclear  spin  selection 
rules  (see  Table  I).  Tf  the  barrier  to  rotation  cin  ages  upon  Sj  -  SQ  excita¬ 
tion,  the  energies  of  the  2e9g  -  2e2g  and  4e2g  -  4e2>,  transitions  would  be 
different;  if  the  magnitude  of  the  change  were  small,  a  doublet  would  be 
produced  in  the  spectrum.  Thus  an  intense  doublet  in  the  spectrum  most  prob¬ 
ably  corresponds  to  the  allowed  2e9g  -  2e2g  and  4e9g  ♦-  4e2g.  transitions. 

A  reasonable  fit  to  the  0°  spectrum  is  obtained  by  assuming  that  in  the 
cluster  ground  state  the  N2  molecule  can  freely  rotate  around  the  benzene-\'2 
bond  and  in  the  excited  state  the  rotation  is  slightly  hindered  by  a  ca.  20 
cm-1  potential  barrier.  The  final  calculated  spectrum  shown  in  Fig  3a  is  a 
superposition  of  both  one-dimensional  rotor  and  vdW  vibrational  transitions: 
due  to  the  rotational  degree  of  freedom,  the  t_,  torsion  is  not  taken  into 
account,  in  the  vibrational  analysis.  Specific  assignments  of  the  calculated 
and  experimental  spectra  are  given  in  Table  II.  As  can  be  seen  in  Fig.  3a  and 


Table  II  the  agreement  between  Hie  experimental  and  calculated  spectrum  is 
very  good  The  fit  is  considerably  worsened  if  the  rotational  barrier  in 
either  state  is  changed  more  than  10  cm  ^ 

T‘  f>(*  spectrum  is  easier-  to  analyze  since  in  this  case  the  origin 
transition  is  allowed  (compare  Table  I)  The  strongest  feature  in  the  sp<- 
tram  in  Fig.  3b.  shifted  bv  ('<  cm  *  with  respect  to  the  benzene  OH  transition, 
is  unambiguously  ascribed  to  the  cluster  F>^  origin.  A  single;  origin  is  ohser  - 
ved  in  agreement  with  the  single  calculated  cluster  geometry  shown  in  Fig.  I. 
The  spectral  features  observed  to  the  red  of  the  origin  are  assigned  as 
hotbands  since  their  intensity  changes  considerably  with  the  carrier  gas 
backing  pressure  but  not  with  ionization  energy.  Several  relatively  weak 
peaks  observed  to  the  blue  of  the  cluster  6*  origin  are  vdW  vibrational  and 
rotational  bands:  their  assignments  are  given  in  Table  II.  The  6*  spectrum 
is  more  than  an  order  of  magnitude  more  intense  tha.  the  0°  spectrum,  with 
relatively  intense  vdW  vibrations  built  on  the  allowed  6*  origin.  The  rota¬ 
tional  transitions  which  are  seen  in  the  0°  spectrum  are  almost  entirely 
hidden  under  the  6^  origin  spectrum.  Moreover,  in  the  region  ca.  lo  cm  1  to 
the  blue  of  the  cluster  origin  (around  the  benzene  6f*  transition),  where  some 
of  the  rotational  transitions  are  calculated  to  appear,  the  mass  detector  is 
still  saturated  by  the  ionized  benzene  molecule  signal.  This  effect  elimin¬ 
ates  observation  of  any  weak  features  and  is  visualized  by  a  dip  in  the 
spectrum.  Further  to  the  blue  only  weak  rotational  transitions  may  be  present 
along  with  vdW  vibrations  as  indicated  in  the  6^  spectrum  in  Fig.  3b. 

Calculations  of  the  rotational  barrier,  rotational  energy  levels,  ami 
vdW  vibrations  confirm  that,  the  benzenefNg^  cluster  0°  spectrum  must  be 
assigned  as  a  superposition  of  both  one-dimensional  rotor  and  vdW  vih: mtiona! 
transitions.  This  approach  gives  the  best  and  most  satisfactory  assignment  of 


the  spectral  features  in  agreement  with  the  selection  rules  given  in  Table  [. 
In  order  to  reproduce  more  accurately  the  most  intense  features  in  the  spec¬ 
trum  a  small,  ca.  20  <m~‘.  barrier  to  rotation  in  the  excited  electronic  state 
is  assumed 

The  duster  ene.gi  shifts  is  ompared  to  transitions  in  benzene  arc  an 
cm-'  and  -f>  cm  1  for  the  r._|  and  ■[.(•  R  transitions,  respectively  This 
unexpected  difference  in  shifts  must  arise  from  local  changes  of  polarizabil¬ 
ity  upon  excitation  of  tin-  benzene  and  can  only  partially  be  explained  by  the 
change  of  the  benzene  fi‘  vibration  in  tin:  cluster.  At.  present  we  do  not  have 
a  satisfactory  explanation  for  these  apparently  different  cluster  shifts. 

B.  Benzene ( CO2) j - 

The  0°  and  (V  spectra  of  the  benzene(C02) j  cluster  are  shown  in 
Fig  4.  The  cluster  energy  shifts  are  -I  cm-1  and  4  cm-1,  for  the  0°  and  6^ 
transitions,  respectively.  One  cluster  configurati  a,  almost  identical  to  the 
nitrogen  case  (shown  in  Fig  1),  is  calculated  for  '‘enzene(C09) j  with  the 
binding  energy  -868  cm-1. 

The  arguments  given  above  for  the  assignment  of  benzene(N2 ) j  cluster 
are  valid  in  the  present  case,  since  both  cluster  symmetries  and  selection 
rules  are  identical.  The  specific  assignments  of  both  spectra  in  Fig.  4  are 
given  in  Table  III.  The  C02  barrier  for  rotation  in  the  excited  state  is 
taken  to  be  16  cm  *  and  good  fit  to  the  0^  experimental  spectrum  is  obtained 
as  shown  in  Fig.  4a.  The  relatively  smaller  energies  of  particular  rotational 
transitions  of  the  benzene(C02)  ^  cluster  compared  to  the  benzene(N’2) .  cluster 
arise  because  the  moment  of  inertia  for  C02  (B  =  1.547  cm-1)  is  larger  than 
for  N ? . 

The  6^  spectrum  in  Fig.  4b  shows  a  strong  origin  peak  with  several  weak 
transitions  to  the  blue.  These  weak  transitions  can  again  be  assigned  (Table 


Ill)  as  internal  rotations  and  vdW  vibrations.  The  origin  corresponds  to  a 
single  calculated  cluster  geometry. 

C.  Benzene (CO) j . 

A  single  cluster-  geometry  shown  in  Fig,  5  is  calculated  for  the 
benzene(CO),  cluster  with  the  binding  energy  -612  cm  1  The  TO  molecule  is 
located  3.24  A  above  the  benzene  ring  forming  a  15"  angle  with  the  benzene 
molecular  plane  in  such  a  way  that  the  oxygen  lies  closer  to  benzene  than  does 
the  carbon  atom.  The  rigid  cluster  symmetry  is  Cs  instead  of  the  C2v  symmetry 
found  for  N'0  and  CO2  clusters  The  CO  rotational  axis  is  no  longer  parallel 
to  the  benzene  six-fold  axis,  the  axes  forming  a  15'  angle  with  respect  to  one 
another.  Therefore,  if  an  internal  rotation  is  present,  the  CO  rotational 
axis  must  be  undergoing  a  small  processional  motion  in  order  to  avoid  a  high 
potential  barrier  to  rotation  The  appropriate  selection  rules  for  the  rota- 
tional  transitions  are  based  ori  the  Cj  2  molecular  s  mmetry  group  (see  Table 
I)  • 

Assignment  of  the  benzcne(CO)j  spectra  is  similar  to  the  assignment  of 
the  benzene(N2)j  spectra,  except  for  the  above  mentioned  selection  rules.  The 
calculated  and  experimental  spectra  are  compared  in  Fig.  6,  with  their  speci¬ 
fic  assignments  given  in  Table  IV.  The  benzene(CO)^  spectra  are  similar  to 
the  analogous  benzene(N2)j  spectra  since  both  these  clusters  are  of  identical 
masses,  have  very  similar  geometries,  and  N2  and  CO  have  almost  Identical 
rotational  constants  (B  =  1.915  cm-1  for  CO  and  1.917  cm-1  for  N2).  Conse¬ 
quently,  both  clusters  exhibit  similar  vdW  vibrations  and  have  similar  ener¬ 
gies  fur  rotational  transitions;  however,  fewer  and  broader  bands  are  observed 
in  the  benzene(CO) ^  spectrum  in  Fig.  6a  compared  to  the  spectrum  of  benzene(S'2l 
in  Fig.  3a,  due  to  an  increased  number  of  allowed  transitions  for  benzene(fO)j 
clusters  (Table  I).  The  overlap  of  transitions  is  also  responsible  for  the 


change  of  the  relative  intensities  in  the  Qo  benzene(ro) ,  spectrum  as  compared 
to  the  0°  spectra  of  the  analogous  clusters  of  S9  am!  C09. 

The  61  spectrum  in  Fig.  fib  is  well  assigned  as  a  superposition  of 
rotational  and  vdW  vibrational  mode  transitions  (Table  IV).  The  spectrum  is 
dominated  by  a  strong  origin  and  vdW  vibrations,  although  weak  rotational 
transitions  can  also  be  .1  i  st  i  ngu  i  shed. 

V.  BENZENE (SOLVENT) 2  CLUSTERS. 

The  analysis  of  the  benzene(solvent)2  cluster  spectra  differs  from  the 
detailed  approach  taken  for  the  benzene(sol vent )|  cluster  spectra  presented 
above  since  the  benzenefsolvent)?  spectra  are  less  well  resolved  and  thus  do 
not  contain  sufficient  experimental  information  about  possible  vibrations 
and/or  rotations.  The  analysis  of  the  benzene(solvent)2  cluster  spectra  is 
based  primarily  on  calculations  of  cluster  geometries  and  assignment  rules 
determined  in  our  previous  studies1”6  of  other  vdW  lusters. 

A.  Benzene (N2 )2 • 

Two  ground  state  minimum  energy  cluster  geometries  are  calculated 
for  the  benzene(N9)?  cluster  (Fig.  7).  The  isotropic  configuration  with 
binding  energy  of  -1007  cm-1  features  a  nitrogen  molecule  on  either  side  of 
the  aromatic  ring.  In  the  anisotropic  geometry  both  nitrogen  molecules  are 
attached  to  one  side  of  the  benzene  ring  with  the  binding  energy  -962  cm”1. 

No  other  local  minima  are  found  in  the  potential  energy  calculations. 

The  benzene(N2)2  cluster  spectra  are  not  detected  in  the  0°  region. 

This  symmetry  induced  transition  (forbidden  in  benzene)  is  certainly  very  weak 
for  the  isotropic  cluster  geometry,  since  the  six-fold  symmetry  of  benzene  may 
be  preserved  in  the  cluster  in  some  averaged  sense,  especially  considering  the 
rotations  of  the  N2  molecules.  Why  the  anisotropic:  cluster  spectrum  is  not 
observed  is  difficult  to  say.  The  61  spectra  taken  at  two  different  ioniza- 


tion  energies  .tre  shown  in  Fig  R.  One  relatively  broad  peak  shifted  -16  cm"1 
to  the  red  of  benzene  transition  dominates  both  spectra  Another  peak  at  9 
cm  *  visible  as  a  shoulder  on  the  main  feature  increases  in  intensity  at  lower 
ionization  energies  ~h  i  >.  feature  does  not  correspond  to  features  in  other 
mass  channels.  The  two  foliates,  one  at  -16  cm  *  and  the  other  at  -9  cm  *. 
are  thus  assigned  to  two.  different  cluster  geometries.  The  cluster  whose  6* 
transition  is  at  -16  tm’1  apparent iy  has  a  higher  ionization  energy. 

i  p 

Benzene  vdW  clusters  with  other  solvents  '  show  that  the  energy  shift 
for  the  isotropic  cluster  is  larger  than  the  energy  shift  for  the  analogous 
anisotropic  cluster.  This  is  due  to  a  better  overlap  of  the  solvent  molecules 
with  the  aromatic  x  cloud  of  benzene  in  the  isotropic  cluster,  resulting  in 
larger  polarizability  changes  upon  cluster  excitation.  Typically,  the  aniso¬ 
tropic  benzenefsol vent )?  cluster  exhibits  a  shift  roughly  equal  to  the  one 
found  for  the  henzeneisolvent)^  cluster.  Given  t h i --  assumption,  the  more 
shifted  peak  in  the  spectrum  is  assigned  to  the  symmetrical  (isotropic) 
cluster  and  the  less  shifted  peak  to  the  anisotropic  cluster.  Comparing  the 
relative  intensities  of  both  features  in  the  spectrum  one  concludes  that  the 
concentration  of  the  benzene!^)?  isotropic  cluster  in  the  beam  is  greater 
than  of  the  anisotropic  cluster  by  roughly  a  factor  of  2  or  3. 

B.  Benzene (COg )2 ■ 

The  benzenefCOg^  cluster  spectra  are  detected  in  both  the  0°  and 
regions  and  are  shown  in  Fig.  9.  The  origin,  shifted  in  both  spectra  -22 
cm  1  to  the  red  of  the  benzene  transition  (at  0  cm"*),  is  ascribed  to  one 
cluster  geometry.  Clearly,  another  cluster  geometry  gives  rise  to  a  peak  at 
-2  cm  5  iri  the  o”  spectrum  in  Fig.  9a.  This  feature  (geometry)  is  not  seen  in 
the  6^  spectrum  due  to  the  saturation  of  the  mass  detector  by  ionized  benzene 
precisely  in  this  spectral  region  (see  dip  in  the  spectrum  in  Figure  9b).  The 


posit  ion  of  t.hp  second  feature  at  -0  cm  1  in  the  dj  spectrum  corresponds  to 
the  origin  of  the  benzenelCO.,)^  cluster-  (compare  with  Figs  11  and  12)  nnd 
thus  is  ascribed  t.o  dissociation  of  rhis  <  luster  giving  'ise  to  i  signal  in 
the  benzene  { CO,,),  mass  ha.nnel.  Pronounced  change  of  the  relative  intensity 
of  this  bar..!  with  varying  the  ionization  energy  confirms  this  assignment 

In  agreement  with  the  experimental  observation,  t-wo  ground  state  mini 
mum  energy  duster  geometries,  simi’ar  to  those  found  for  benzene) S', )  ,  c. luster 
(shown  in  Fig.  7)  are  >•  c  1 ;  a;  1  a  I  ed  The  binding  energies  are  -1716  < .m  1  and 
-1955  cm"1  for  the  isotropic  and  anisotropic  clusters,  respectively,  The  -22 
cm-1  shifted  origin  in  the  two  spectra  is  assigned  to  the  isotropic  cluster, 
whereas  the  2  cut”'  shifted  origin  in  the  0°  spectrum  is  assigned  to  the 
anisotropic  cluster.  As  can  be  seen  front  the  comparison  of  the  relative 
intensities  in  Fig.  Da.  the  concentration  of  the  isotropic  cluster  in  the 
supersonic  molecular  jet  is  larger  than  the  concen t r i t i on  of  the  anisotropic 
cluster  The  intensity  of  the  peak  ascribed  to  the  misotropic  cluster  can  be 
diminished  by  increasing  the  ionization  energy  Tl.  behavior  of  the  ben- 
zene(C02)2  cluster  is  thus  very  similar  to  the  one  reported  above  for  the 
benzene (.V,,) £  cluster. 

C.  Benzene (CO) 2* 

Four  cluster  geometries  (one  isotropic  and  three  anisotropic)  arc 

calculated  for  the  benzene(C0)2  cluster  with  the  binding  energies  -1284  cm"1 

and  -1580  cm  1 ,  -1601  cm1  and  -1697  cm-1,  respectively. 

The  61  spectrum  of  the  benzene(CO),  cluster  (Fig  10)  is,  different 

from  the  spectra  of  analogous  clusters  with  \r2  and  CO ,.  The  spectrum  is 

1  -1  _ 

relatively  weak  and  exhibits  four  distinct  features  at  -56  cm”1.  -61  cm  1 

cm  1  and  -88  cm-1  with  respect  to  the  benzene  origin.  The  strongest  feature 
at  -56  cm"1  probably  corresponds  to  the  isotropic  calculated  cluster  geometry 


which  is  found  to  occur  with  the  largest  probability.1'  TllH  other-  features 
may  be  assigned  to  the  other  three  > a  1 cu lat ed  anisotropic  cluster  geometries. 

The  broad  background  in  the  spectrum  is  likely  due  to  dissociation  of 
higher  order  -'lusters  hut  cannot  tie  entirely  eliminated  by  lowering  the  ion 
i/.ation  energy.  This  may  indicate  that  some  of  the  peaks  in  the  spectrum  are 
due  t.o  dissociation  of  -  lust  -rs  instead  of  being  inherently  related  to 

the  benzeoetCO),  clusters  Higher  order  clusters  of  benzene  with  CO  re  not. 

however,  detected  in  *  hose  ■•xper  i  men  t  s . 

VI.  LARGE  CLUSTERS  -  TRANSITION  ENERGY  SHIFTS. 

Except  for  the  hen/ene ; C09 cluster  spectrum  shown  in  Fig.  11,  spectra 
of  large  clusters  are  not  detected  in  the  0^  region.  The  benzene  0°  transi¬ 
tion  is  cluster  symmetry  induced  and  thus  is  not  expected  to  be  intense  for 
low  concentrations  of  specific  clusters.  One  origin  shifted  by  -9  cm  1  is 
observed  in  the  0*’  spectrum  of  the  benzene (COg S3  cl  st.er. 

The  6 1  spectra  of  benzene  solvated  by  up  to  gbt  N’2  molecules  and 
seven  C0o  molecules  are  shown  in  Figs.  12  and  13,  respectively.  Progressive 
transition  energy  shifts  toward  lower  energies  with  increasing  cluster  size 
are  observed  for  the  large  clusters.  One  rather  broad  peak  dominates  the 
spectra  of  N?  clusters,  whereas  spectra  of  progressively  larger  C02  clusters 
are  also  generally  broad  but  do  show  some  structure  on  a  broad  background. 

The  C02  cluster  spectra  are  in  addition  somewhat  obscured  by  a  pronounced  dip 
due  to  the  saturation  of  the  mass  detector  by  ionized  nonclustered  benzene 
molecules.  Plots  of  an  average  cluster  transition  energy  shift  as  a  function 
of  cluster  size  are  given  in  Fig.  14.  A  linear  dependence  of  the  energy  shift 
on  the  cluster  size  is  observed  for  nv2.  Addition  of  each  solvent  molecule 
thus  causes  a  similar  perturbation  of  the  benzene  transition  once  the  first 
two  solvent  molecules  are  attached.  For  clusters  with  seven  and  eight  nifro- 


gen  molecules,  the  change  in  transition  energy  shift  seems  to  be  slightly 
smaller,  possibly  due  to  a  gradually  weaker  interaction  of  the  solvent  mole¬ 
cules  with  the  aromatic  ■  cloud  of  benzene  as  the  duster  size  increases 

The  above  results  suggest  that  the  cluster  shifts  are  not  entirely 
saturated  by  one  solvent  molecule  on  each  side  of  the  aromatic  ring,  as  is 

t  P 

found  for  benzene-alkane'  and  t u 1 u»-ne -a  1 kane  clusters.  Larger  shifts  still 
arise  with  addition  of  more  than  two  molecules  indicating  that  several  solvent 
molecules  may  effectively  interact  with  benzene.  This  is  certainly  the  case 
for  the  solute-solvent  systems  studied  here  for  which  the  sizes  of  linear 
solvent  molecules  are  relatively  small  compared  to  benzene.  The  N'g  clusters 
with  benzene  exhibit  larger  energy  shifts  than  the  corresponding  clusters  of 
fC2,  even  though  the  CO?  molecule  has  a  relatively  larger  polarizability.  The 
increase  of  the  energy  shift  with  cluster  size  is  also  larger  for  benzene(N'2) 
than  for  benzene(C02)r|.  This  is  probably  due  to  the  smaller  size  of  the  N2 
molecule  and  thus  a  better  spatial  packing  of  nitn  .<-n  around  the  benzene  and 
concomitant  stronger  interaction  with  the  electric  •  cloud. 

Many  solvent  molecules  effectively  interact  with  benzene  and  contribute 
to  the  cluster  transition  energy  shift  in  large  clusters.  Since  only  one 
rather  broad  feature  is  observed  in  the  large  cluster  spectra,  the  spectral 
energy  differences  between  various  cluster  configurations  must  he  relatively 
small.  Comparison  of  the  spectra  of  benzene(»2)n  clusters  with  the  spectra  of 
benzene(C02)n  clusters  suggests  that  the  N2  clusters  exhibit,  fewer  and  more 
unique  orientations  than  the  C02  clusters;spectra  of  N’2  clusters  show  one 
relatively  narrow  peak,  while  those  of  C02  are  broader  and  have  some  structure 
on  the  main  background.  This  difference  may  be  due  to  the  larger  size  of  the 
C02  as  compared  to  \’2.  Calculations  of  the  possible  large  cluster  geometries 
are  presently  being  carried  out  in  order  to  confirm  this  interpretation 
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Preliminary  results  indicate  that  clusters  with  more  than  two  solvent  mole¬ 
cules  exhibit  a  large  number  of  local  potential  energy  minima  differing 
slightly  in  cluster  binding  energy.  The  number  of  local  minima  increases  with 
the  cl  ust.er  size 

If  large  clusters  an*  to  serve  as  potential  models  of  condensed  phase 
systems,  gas  to  cluster  energy  shifts  must  be  explored  and  compared  with  those 
of  condensed  phase  solutions.  The  largest  energy  shift  observed  in  our  exper¬ 
iments  for  benzeneiMg),,  clusters  is  ca.  -42  cm"1  for  a  cluster  with  eight  N’p 

1  IQ 

molecules.  A  shift  of  about  50  cm  is  measured  for  benzene  dissolved  in 
supercritical  nitrogen  fluid  (in  the  gas  phase)  at  a  fluid  density  of  ca.  220 
kg  m  3  (at  T  =  205  K  and  P  -  ca.  300  bar).  This  shift  increases  with  increas¬ 
ing  fluid  density  (pressure)  due  to  forced  crowding  of  the  solvent  around  the 
solute  so  that  more  molecules  effectively  interact  with  benzene.  In  the 
liquid  phase  and  high  density  fluid  phase  the  energy  shift  increases  up  to  - 
180  cm-1  due  to  the  additional  interaction  of  nitre  ,en  molecules  forced  to 
crowd  about  the  benzene  molecule  at  its  periphery 

Based  on  the  comparison  of  benzene(N2)n  cluster  shifts  with  the  super¬ 
critical  nitrogen  fluid  results,  one  concludes  that  in  the  low  density  conden¬ 
sed  phase  "solvation"  of  benzene  takes  place  mostly  about  the  ,t  system  and  not 
at  the  more  repulsive  C-H  periphery  positions.  This  latter  type  of  solvation 
must  be  forced,  however,  in  high  density  fluids  and  the  liquid  state. 

VII.  NUCLEATI0N. 

The  nucleation  process  for  small  vdW  clusters  can  be  largely  understood 
only  if  various  spectral  features  can  be  ascribed  to  specific  cluster  geo¬ 
metries.  Two  basic  types  of  nucleation  can  be  distinguished:  homogeneous 
nucleation  in  which  solvent  molecules  are  added  to  the  solute  molecule  or 
cluster  one  molecular  at  a  time;  and  inhomogeneous  nucleation  in  which  more 
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than  one  solvent  molecule  is  added  to  the  solute  or  cluster  at  a  given  time. 
Relative  intensity  data  for  clusters  of  benzene  and  toluene  with  hydrocarbon 
solvents  has  led  to  the  conclusion  that  solvent  molecules  exist  in  the  super¬ 
sonic  expansion  in  the  form  of  dimers  or  larger  aggregates  and  that  isotropic 
benzene( so  1  vent  ) 2  clusters  ale  homogeneously  nucleated  and  anisotropic 

1  O 

clusters  are  i nhomogeneous 1 y  nucleated.  Whether  homogeneous  nucleation 

(isotropic  clusters)  or  inhomogeneous  nucleation  (anisotropic  clusters)  arises 
deperds  on  the  relative  size  of  the  solvent  dimer  binding  energy  ( E , )  with 
respect  to  the  solute  solvent  (cluster)  binding  energy  (F.  )  in  the  cluster 
The  smaller  the  binding  energy  of  a  solvent  dimer  compared  to  the  binding 
energy  of  the  cluster,  the  higher  the  concentration  of  isotropic  clusters. 
Inhomogeneous  nucleation  is  found  to  predominate  for  small  clusters  of  benzene 
and  toluene  with  two  hydrocarbon  solvent  molecules  because  E^  Ec  is  relatively 
large  (0.3  -  0.8  depending  on  the  solvent). 

The  binding  energies  of  N2  and  C02  dimers  at--,  however,  quite  small: 
-150  cm-1  and  -442  cm'1,  respectively.  As  shown  in  Table  V  the  Er  value 
for  the  benzene(N'2)2  and  benzenelCC^^  clusters  is  only  0.15  and  0  23,  respec¬ 
tively.  This  clearly  favors  homogeneous  nucleation  and  formation  of  isotropic 
clusters,  since  solvent  dimers  can  easily  be  dissociated  with  very  little 
excess  energy  Our  spectra  confirm  this  Implication  because  in  all  cases 
studied  the  peak  ascribed  to  isotropic  benzene(solvent)2  clusters  is  much  more 
intense  than  the  peak  due  to  anisotropic  clusters.  This  effect  is  especially 
strong  in  the  case  of  N’2  clusters  for  which  Ecj/Ec  is  small. 

Clusters  larger  than  benzene(solvent)2  are  most  likely  formed  by 
inhomogeneous  nucleation.  First,  larger  solvent  clusters  have  larger  binding 
energy.  High  solvent  aggregate  binding  energy  favors  inhomogeneous  mnlea- 
tion.  Comparison  of  the  calculated  solvent  dimer  and  trimer  binding  energies 


given  iri  Table  V  shows  an  almost  threefold  increase  of  the  solvent  aggregate 
binding  energy  ( E a )  for  the  it  inner  with  respect  to  the  dimer.  At  the  same 
time  only  a  twofold  increase  of  the  cluster  binding  energy  takes  place  for  the 
benzene  (solvent.^  cluster  compared  to  the  benzene  (solvent  )j  cluster  which 
results  in  an  increase  of  the  E  i\  ratio  with  the  increase  of  cluster  size 
This  strongly  favors  i  nhomogeneoc  -  iu< .  I  eat.  ion  for  large  benzene(so  I  vent. ) 
clusters.  Second,  as  the  'luster  ^rows  in  size  the  number  of  bonds  over  which 
to  distribute  the  binding  energy  of  r he  collision  partner  increases  greatly, 
and  thus  the  ability  of  the  cluster  to  dissipate  the  added  binding  energy 
i ncreases . 

VIII.  CONCLUSIONS. 

Two  color  time  of  flight  mass  spectroscopy  has  been  employed  to  study 
small  and  large  clusters  of  benzene  with  nitrogen,  carbon  monoxide,  and  carbon 
dioxide  created  in  a  supersonic  molecular  jet.  P o'  utial  energy  calculations, 
normal  coordinate  and  internal  rotational  analyses  h ive  been  employed  for  the 
assignment  of  the  0®  and  spectra  of  the  benzene! so  1  vent ) j  clusters.  Geo 
metries  of  small  clusters  (n=l,2)  have  been  computer  calculated  and  assigned 
to  experimental  spectra.  Transition  energy  shifts  for-  both  small  and  large 
clusters  with  up  to  eight  solvent  molecules  are  investigated  for  the  first 
time  as  a  function  of  cluster  size.  The  following  major  conclusions  emerge 
from  this  study: 

1.  The  calculated  cluster  binding  energies  scale  well  with  the  solvent 
polarizabilities;  the  larger  the  polarizability  the  higher  the  binding  energy. 

2.  The  cluster  spectral  shifts  depend  on  the  solute-solvent  geometry 
•ad  interaction  but  are  less  dependent  on  solvent  polarizability.  The  depend¬ 
ence  of  the  energy  shift,  on  the  cluster  geometry  is  found  to  he  less  prominent 
than  in  our  previc us  studies  of  benzene-alkane  clusters.  Both  red  and  blue 


cluster  energy  shifts  are  observed  and  the  additive  energy  shift  rule  for  the 
isotropic  benzene-solvent  clusters  does  not  apply  in  general  for  the  clusters 
stud  ied . 

3.  Analysis  of  the  rotational  and  vibrational  motions  within  the 
benzenef so i ven t )  |  clusters  makes  possible  the  assignment  of  the  rotational  and 
vdW  vibrational  modes  in  the  spei I . a  Both  internal  rotation  modes  of  the 
solvent  molecule  rotating  mound  the  solute-solvent  bond  axis  and  vdW  vibra¬ 
tions  ai e  active  in  the  benzene! so  1 ven t ) j  cluster  spectra  and  must  be 
accounted  for  All  three  solvent  molecules  studied  rotate  nearly  freely 
around  the  benzene- so  1 ven l  bond  axis  The  spectra  are  well  assigned  assuming 
no  barrier  to  rotation  in  the  ground  state  (free  rotor)  and  a  small  ca.  2D 

cm  !  harrier  to  rotation  in  the  ex>  ted  state  (slightly  hindered  rotor).  The 
solvent  rotational  axis  for  these  linear  diatomic  and  ti. atomic  molecular 
solvents  is  oriented  along  the  benzene-solvent  bond  and  thus  little  or  no 
barrier  for  the  internal  rotation  exists. 

4.  The  benzenefso 1  vent ) j  vdW  vibrations  obs  i ved  are  those  involving 
bending  motions  parallel  to  the  benzene  molecular  plane  (bx  and  by),  a 
stretching  motion  along  the  solute-solvent,  bond  (Sz),  and  a  torsional  motion 
(ty).  The  t?  motion  is  not  present  because  of  the  free  rotation  around  the  z 
axis 

5.  Differences  pertaining  to  energy  shifts,  intensity  distribution  and 
appearance  of  the  rotational  and  vdW  modes  are  found  between  the  O1^  and 
cluster  spectra.  These  differences  arise  mostly  from  different  selection 
rules:  the  origin  is  allowed  at  the  6^  transition  but  forbidden  at  the  0®. 
Consequently,  relatively  weak  internal  rotational  mode  structure  is  identified 
in  the  0”  spectra  but  is  difficult  to  distinguish  in  the  6*  spectra  which  are 
dominated  by  a  strong  origin  and  associated  vdW  vibrations 


(>.  Large  vdW  i  lusters  nf  up  to  eight  molecules  and  seven  C0-, 
molecules  solvating  benzene  are  observed  for  the  first  time  One  rather  broad 
cluster  feature  dominates  the  spectra.  A  I i near  energy  shift  to  the  red  with 
increasing  cluster  si/e  is  found  The  cluster  shifts  are  not  saturated  by  one 
solvent  molecule  on  each  side  of  the  aromatic  ring  as  was  found  previous.  for 
small  benzene -,t  1 kane  clusters  This  may  be  due  to  the  relatively  small  size 
of  the  solvent  molecules. 

7.  Homogeneous  nuc  leaf  ion  is  found  to  dominate  the  formation  of  small 
benzenet so  1 ven t )2  clusters.  This  is  due  to  a  very  small  solvent  dimer  binding 
energy  with  respect  to  the  cluster  binding  energy.  In  large  clusters  the 
ratio  of  the  solvent  aggregate  binding  energy  to  the  cluster  binding  energy 
and  the  number  of  ways  to  share  the  binding  energy  between  the  cluster  bonds 
increase  considerably;  hence,  inhomogeneous  nucleation  may  be  favored  in  the 
formation  of  large  clusters. 
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TABLE  I 


Rovibronic  and  nuclear  spin  selection  rules. 

*  and  -  stand  for  an  allowed  and  forbidden  transition,  respectively. 


Transitions  not 

indicated  are 

completely 

forbidden . 

C 1  uster 

Molecular 

Symmetry 

Group 

Possible 

Transitions 

Selection 

Rovibronic 

0°  6 1 

0  0 

Rules 

Nuclear  Spin 

"ig  - 

"  alg 

- 

4- 

* 

benzene (N2 ) j 

e2g  ~  e2g 

4- 

4- 

- 

& 

G24 

alg  “ 

~  P2g 

4- 

* 

- 

benzene(C02)2 

plu  ~  elu 

4- 

4- 

- 

b2u  ~  b2u 

- 

- 

- 

elu  " 

"  b2u 

- 

- 

- 

al  - 

-  al 

- 

* 

4- 

e2  ^ 

-  e2 

4- 

4- 

♦ 

b  ne(C0)j 

G12 

al  ~  e2 

4- 

- 

4- 

el  - 

^  el 

4- 

- 

4- 

b2  *— •  b2 
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TABLE  II 


Assignments  of  the  0^  and  6*  benzene  ( \’2  )  j  cluster  spect  ra 
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Transition 

Region 


TABLE  iff 


Assignments  of  the  0^  and  6^  benzene ( CC^ ) j 
r  luster  spectra  shown  in  Fig.  4. 


Exper imenta 1 
Peak  Positions 
(Relative  to 
cluster  origin! 


Calculated 
Rotational 
Trans i t ions 


Ca  1  <u  1  a  t  ed 
vdW  Vi  bra  t ions 


-  1 
cm  * 

Value 
cm- 1 

Ass  igriment 

Value 
cm' 1 

-8 

-7 

2p2g  “ 

0alg 

- 

-1 

-1 

2e2g  - 

2e2g 

- 

o  c 
c 

1 

1 

4e2g  *• 

4e2g 

- 

(Fig.  4a) 

7 

6 

0alg  - 

2e2g 

- 

[38.085  cm'1  ] 

17 

25 

19 

26 

2e2g  ^ 

0alg  *■ 

4p2g 

;,>2g 

- 

33 

- 

35 ;  38 

65 

_ 

72;  70 

Ass i gnment 


o  9 

V;V 


sz:bx 
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TABLE  TV 


Assignments  of  the  0^  and  61  benzene(C0Jj 
cluster  spectra  shown  in  Fig.  6. 


Trans i t ion 

Experimental 

Calculated 

Calcu  1  at.ed 

Region 

Peak  Positions 

Rotat ional 

vdW  Vibrations 

(Relative  to 

Trans i t i ons 

cluster  origin) 

Value  Assignment 

Value  Assignment 

cm” 1 

cm'  1 

cm- 1 

-23 

-24 

4e2 

-  2e2 

- 

- 1 5 

-16 

3b2 

-  Oa^ 

- 

-8 

-8 

2e2 

-  Oa^ 

- 

o  o 
o 

0 

-1 

2e2 

-  2e2 

- 

(Fig.  6a) 

0 

0.5 

3b2 

-  3b2 

- 

[38.110  cm'1 

0 

1 

4e0 

-  le2 

_ 

10 

7 

Oaj 

-  2e2 

10 

1 1 

Oaj 

-  3b2 

14 

bx 

23 

21 

Oa^ 

-  3b2 

18 

by 

23 

24.5 

2e2 

-  4e2 

35 

33 

0a1 

-  4e2 

42 

- 

- 

42 

bx3 

55 

54 

ty.b/ 

66 

68 
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TABLE  IV  -  Continued  .2  of  2. 


Transition 

Experimental 

Calculated 

Region 

Peak  Positions 

Rotational 

(Relative  to 

Trans i t ions 

cluster  origin) 

Value  Assignment 

cm-1 

cm'  * 

Ca leu  1  a  ted 
vdW  Vibrations 


Value  Assignment 

cm- 1 


Figure  2 


Figure  3 


FIGURE  CAPTIONS 


Figure  1  Calculated  ground  state  minimum  energy  configuration  (a)  and 

eigenvalues  and  eigenvectors  of  the  vdW  modes  (b-f)  for  hen- 
zene(NT2>|.  Rigid  cluster  symmetry  is  taken  to  be  C2v 


Energies  of  the  internal  rotational  levels  of  the  N2  molecule  in 
the  benzenefN, ! ^  cluster  calculated  with  B  =  1.917  cm-1  as  a 
function  of  a  Vfi  potential  barrier.  The  symmetries  of  the  free 
rotor  rotational  levels  at  0  cm-1,  in  order  of  increasing  energy, 
are  0alg,  lelu,  2e2g,  3b2u.  4e2g,  5eU),  6alg  and  7e]u. 


Two  color  TOFMS  of  the  benze»e(N2)j  (luster  in  the  0°  (a)  and 
6o  (b)  regions.  The  internal  rotational  (continuous  line)  and 
vdW  vibrational  (dashed  line)  transitions  are  calculated  as 
described  in  the  text  with  the  potential  barrier  in  the  excited 
state  Vg  (Sj)  =  20  cm-1.  The  0  energy  in  the  0°  spectrum  corres¬ 
ponds  to  the^  calculated  position  of  the  forbidden  origin.  The 
intensities  in  the  calculated  spectra  are  chosen  arbitrarily  in 
agreement  with  the  experimental  spectra.  Specific  assignments 
are  given  in  Table  II.  The  arrows  indicate  positions  of  the  0^ 
(38086.1  cm  *)  and  6*  (38608.5  cm-1)  transitions  in  bare  benzene. 


V 
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Figure  4 


Two  color  TOFMS  of  the  benzene (C02) j  cluster  in  the  0°  (d)  and  61 


Figure  5 


Figure  6 


Figure  7 


(b)  regions.  Calculated  spectra  (Vfi  (Sj)  -  16  cm'1)  consist  of 
rotational  (continuous  line)  and  vdW  vibrational  (dashed  line) 
transitions  The  0  energy  corresponds  to  the  calculated  position 
of  the  forbidden  origin.  Specific  assignments  are;  given  in  Table 
m.  The  arrows  indicate  corresponding  benzene  transitions  (see 
caption  for  Fig  3).  The  dip  to  the  left,  of  the  origin  in  b  is 
due  to  the  saturation  of  the  mass  detector  by  ionized  benzene 
molecules . 


Calculated  ground  state  minimum  energy  configuration  for  the 
benzene(CO)-j  cluster.  The  rigid  cluster  symmetry  is  Cg.  The  vdW 
vibrations  are  very  similar  to  those  shown  in  Fig.  1  for  the 
benzene^)}  cluster  with  the  selection  rules  given  in  section 
III. A.  The  solid  line  indicates  the  CO  rotational  axis. 

Two  color  TOFMS  of  the  benzene(CO) j  cluster  in  the  0®  (a)  and  61 
(b)  regions.  Specific  assignments  of  both  experimental  and  cal¬ 
culated  spectra  are  shown  in  Table  TV.  The  arrow  indicates  the 
benzene  origin  (compare  Fig.  3). 

Calculated  ground  state  minimum  energy  configurations  of 
benzenelNgJg  clusters:  isotropic  (a)  and  anisotropic  (b). 


Figure  8  Two  color  TOFMS  of  benzene(N2)2  clusters  taken  at  two  different 

ionization  energies:  .'17270  cm  '  (upper)  and  36550  cm’1  (lower). 

0  corresponds  to  the  6^  transition  of  benzene  at  38608.5  cm  '. 

The  most  intense  peak  in  the  spectrum  is  assigned  to  the  iso¬ 
tropic  cluster 

Figure  9  The  benzene!  00.;, ),  two  color  TOFMS  at  0®  (a)  and  6^  (h)  transi¬ 
tions.  The  scale  is  relative  to  the  benzene  transitions. 

Figure  10  Two  color  TOFMS  of  the  benzene(C0)2  cluster  at  The  scale  is 

relative  to  benzene  6^  transition. 

Figure  11  Two  color  TOFMS  of  the  benzene(C02)2  aster  in  the  0^  region. 

Figure  12  The  6j!  two  color  TOFMS  of  the  benzene!. \->)n  clusters.  The  spectra 

are  numbered  according  to  the  number  of  N2  molecules  in  the 
cluster.  0  corresponds  to  the  benzene  6^  transition  at 
38608.5  cm- 1 . 

Figure  13  The  6^  two  color  TOFMS  of  the  berizene(C02)n  clusters.  The 

numbers  reflect  the  number  of  C02  molecules  in  the  cluster.  The 
scale  is  analogous  to  that  of  Fig.  12. 

Figure  14  Clusters  transition  energy  shifts  plotted  as  a  function  of  the 
cluster  size.  The  error  bars  indicate  uncertainty  due  to  the 
broadness  of  the  spectral  features. 
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